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Abstract 

We introduce and study a new class of power-counting non-renormalisable gauge theories in 
four space-time dimensions. The Lagrangian is an arbitrary function of the self-dual part of the 
field strength. The resulting perturbation theory has the property that whenever two derivatives 
act on an internal line propagator, the result is a delta-function and the line collapses to a point. 
This means that there remains at most one derivative on each internal line, which gives improved 
ulta-violet behaviour. For many purposes, this class of theories behaves just like ordinary Yang- 
Mills theory. In particular, they all share the Yang-Mills theory MHV amplitudes. Moreover, these 
theories remain constructible (in the generalised sense), with the higher-point tree level scattering 
amplitudes obtainable from the lower-point amplitudes using the BCFW recursion relations, and 
adding new amplitudes at every particle number. Also, the square of these gauge-theory amplitudes 
gives the scattering amplitudes of “deformation” of General Relativity, at least for the low par¬ 
ticle numbers that we checked. We compute the one-loop beta-function of the first new coupling 
constant, and find it to be positive, which signals the associated non-renormalisable interaction 
becoming important in the ulta-violet. 


1 Introduction 

By now it is commonly appreciated that every field theory is renormalisable provided all terms com¬ 
patible with the symmetries are included into the Lagrangian; see [T] for a recent authoritative account 
of effective field theory. Moreover, it is in practice unnecessary to include all possible terms, as many 
of the terms can be eliminated by field redefinitions that do not change the physics, see e.g. [2] for a 
discussion of this in the context of the renormalisation group flow. 

This (and related [3]) paper can be viewed as an attempt to construct a renormalisable effective 
field theory model where only a very restricted type of terms is allowed to appear in the Lagrangian. 
Thus, we introduce and study a model where the number of terms in the Lagrangian is infinite, but 
much fewer as compared to what would be allowed in the full effective held theory. Our hope is that 
the class of theories we consider is still large enough to be closed under renormalisation. There is a 
simple mechanism, to be described below, that gives a justihcation for our (possibly over-optimistic) 
expectations, but beyond one loop we do not know whether our hope is realised. At one loop a dehnite 
statement can be made; see below. 

Thus, the purpose of this paper is to introduce and study a new inhnite-parametric class of power¬ 
counting non-renormalisable theories. What is new and striking about these theories is that for many 
purposes they behave like the renormalisable theory (mother theory) from which they originated. The 
mother theory here is Yang-Mills theory in four space-time dimensions, and so we will refer to the 
new theories as deformations of Yang-Mills (YM) theory. 
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1.1 The class of theories 


The new theories are obtained by augmenting YM with certain new, importantly chiral, see below, 
power counting non-renormalisable interactions. The zeroth order idea is to consider effective La- 
grangians of the form 
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( 1 ) 


where the first term is the usual YM Lagrangian, and the new terms are gauge-invariant scalars con¬ 
structed from higher powers of the curvature. Note that there can be more than one new coupling 
constant at each order in the curvature, and we wrote the Lagrangian just schematically. Clearly, the 
coefficients in front of the new terms are dimensionful, and it is convenient to represent them as dimen¬ 
sionless couplings a, P,... times an appropriate inverse power of some mass scale M. Alternatively, 
the Lagrangians of interest can be written as 


C = M^f{F/M^), 


( 2 ) 


where / is a gauge- and Lorentz-invariant function of its curvature argument F = where a is the 
Lie algebra index. 

It is important to emphasise that the class of theories that we would like to consider is much 
smaller than the class of general effective field theories for the gauge field. Indeed, the effective field 
theory Lagrangian would also contain terms involving derivatives of the curvature. Thus, another 
gauge-covariant object that can be used in the construction of Lagrangians is the Lie-algebra valued 
vectorY“ := (FF^j^. Building blocks with even larger numbers of derivatives are possible, e.g. d^df^Fp^. 
It is clear that Lagrangians involving powers of Y“ would lead to field equations containing more 
than second derivatives of the gauge field. In contrast, the Lagrangians Q lead to second order in 
derivative field equations. It is clear that some of the higher-derivative terms can be eliminated by 
field redefinitions. But presumably there are terms that cannot be disposed off this way, and so Q is 
too restrictive as an effective field theory Lagrangian. 

Lagrangians of the type Q do appear as effective ones in a variety of theories. For example, one 
can imagine some fermionic fields interacting with the gauge field being integrated out. The resulting 
one loop effective action is non-local, but admits an expansion in powers of the field strength, with 
the scale M being the typical mass of the particles integrated out. The famous example of an effective 
Lagrangian appearing this way is due to Euler-Heisenberg [Ij. Of course, such effective Lagrangians 
typically contain also the derivatives of the field strength, with Q being just the constant field 
approximation. 

Since the Lagrangian Q is non-renormalisable but contains much fewer terms than the effective 
field theory Lagrangian, one should expect that counterterms needed to cancel the arising loop diver¬ 
gences are not necessarily of the type already included in Q. This is indeed the case. In other words, 
having added some non-renormalisable interactions one cannot stop and has to add all of them. The 
result is the effective field theory Lagrangian containing higher derivatives. This is not the route we 
want to follow in this paper. Instead, we will be studying theories of type Q with a certain further 
restriction added. 

We will further cut the class (§ as follows. We have already restricted the theories of interest 
as compared to the general effective field theory by requiring the field equations to be second order 
in derivatives. This requires the Lagrangian to be a function of the field strength. We now further 
restrict what type of functions are allowed by requiring that all of our theories must share the same 
instanton sector. Recall that YM anti-self-dual (ASD) instantons are field configurations that have 
vanishing self-dual (SD) part of the curvature tensor with a being the Lie algebra index. It is 
convenient to describe the self-dual/anti-self-dual decomposition of the curvature in terms of spinor 
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notations. In our conventions 


^MM’NN' — 2^MN^M'N' + 2^M'N'^NM- (3) 

Here Fmm'nn' spinor form of the curvature where each space-time index is replaced by a pair 

of spinor indices, and , F^,j^, = F°'^m'mn' are the self- and anti-self-dual parts of 

the curvature. Instantons, or ASD gauge fields are then field configurations satisfying 

Ft,N = 0- (4) 

It is not hard to see that they are automatically solutions of the Yang-Mills equations. Indeed, as is 
well known, modulo a surface term, the YM Lagrangian can be rewritten as 

Fym = -^{F^Nf ■ (5) 

The field equations following from this Lagrangian are 

dN'^Ft,^ = 0 . ( 6 ) 

It is clear that the ASD gauge fields satisfying Q are also solutions of Q. 

We now require this to continue to be true for an arbitrary member of our family of theories. It is 
easy to see that this implies changing Q into 

£ = M^f{F^j,/M% (7) 

i.e. into a function of only the self-dual part of the field strength. The resulting field equations are 
then 

dN'M ( J!/a ) = 0 - ( 8 ) 

It is clear that on instantons @ the matrix of first derivatives of / appearing in the field equations 
is just a constant (possibly zero), and thus all field configurations Q are also solutions of the field 
equations Q. 

One of the main goals of the present work is to characterise the class of theories 0 in as much 
detail as possible. In particular, we shall soon see that there is a simple mechanism that makes these 
theories behave like the usual Yang-Mills theory, in spite of a clear presence of power-counting non- 
renormalisable interactions. There are many other nice properties that theories Q have, and this is 
what motivates our interest in them. 


1.2 A 2-parameter example 

As a concrete example of a theory of the type 0 let us consider a gauge theory containing just a 
single additional interaction term as compared to the standard YM 
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(9) 


When a perturbative expansion is performed and the coupling constant is absorbed into the connection, 
one obtains ag/M'^ in front of the new term, which is a convenient parameterisation of the new 
interaction, as we shall see. 

An effective field theory model of type Q has been considered in the literature before; see in 
particular [5j for a study in the context of amplitudes. However, for reasons of unitarity, see below, 
the new interaction in Q was always considered with its parity-dual. It is important to stress that 
we only allow in the Lagrangian the chiral half of this F^-type interaction. 
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1.3 A summary of properties 

Here we give a brief account of the properties of 0: with some of these properties derived below, and 
some others to be obtained in [3]. Our main aim here is to indicate how nicely these theories behave, 
and thus hopefully stimulate further interest. 

Our main interest in Q is due to the fact that, in spite of seemingly being a too small class of 
theories to remain closed under renormalisation, there is a simple mechanism that renders the diver¬ 
gences of these non-renormalisable theories essentially those of their renormalisable mother theory. 
To describe this mechanism, let us consider the perturbative expansion of Lagrangians 0. It is clear 
from 0 that the gauge field kinetic term is unchanged from what it is in YM theory. So, our theories 
describe interacting gluons, with some additional chiral interactions added on top of what is present 
in YM. An important observation is that the derivative of the gauge field enters the Lagrangian only 
in a very special combination 

{da)‘XB ■= . ( 10 ) 


Unlike in YM theory, there are now interaction vertices with as many derivatives as the valency of the 
vertex. However, and this is important, there is at most just a single derivative acting on every leg 
of the vertex. So, at first sight our theories appear to diverge much worse than YM. However, there 
is a special property that when two derivatives happen to act on the same internal line propagator, 
the resulting second derivative is a multiple of the box operator, and then everything gets replaced by 
the (5-function, see (55) for a graphical representation. As the result, many internal lines in Feynman 
diagrams of our theory will collapse to form effective higher-valent vertices. After the collapse, only 
internal lines with at most one derivative of the propagator remain. It can be seen that this is also what 
happens in Yang-Mills theory - after some cancellations only diagrams with at most a single derivative 
on an internal line remain. Given that the number of derivatives on internal lines is what determines 
how diagrams diverge, we conclude that the divergences of our class of theories are essentially the 
same as in YM. See below for a further discussion of this, and also for the associated power-counting. 
This is the main justification for our hopes that the class 0> while clearly much smaller than the full 
effective field theory, may still be large enough to be closed under renormalisation. 

This property of our theories is at work at any loop order, but we have not yet carefully analysed 
its consequences beyond the one loop. At one loop order there is a definite statement that can be 
made: The class of theories (j^ is one-loop renormalisable. Thus, in spite of being clearly very far 
from the general effective field theory for the gauge field, no new one loop counter terms that are 
not already contained in (|^ are required. All one loop divergences can be absorbed into either field 
redefinitions or coupling constant renormalisations. There will be plenty of illustrations of why this is 
to be expected in the present paper, but the detailed proof of this statement is spelled out in [3]. 

Another important set of properties concerns the scattering amplitudes for ([^ . It is by now a well- 
known story that Yang-Mills theory is constructible, in the sense that all gluon scattering amplitudes 
can be determined by certain recursion relations from the basic 3-gluon ones. This means that all 
scattering amplitudes can be determined without recourse to any Lagrangian or Feynman rules, simply 
by constructing higher particle number amplitudes from lower ones using the recursion. In this paper 
we shall show that the same statement is true (in a generalised sense) for theories 0. One only has 
constructibility in a generalised sense because in conventional constructible theories like YM or GR 
there is only one coupling constant. As a result, one only needs the 3-point amplitudes as the seed 
for recursion. In our case, there is an infinite number of couplings. As we shall see, at the level of 
recursion relations these new coupling constants will enter at each particle number. Thus, at each 
gluon number a certain set of new amplitudes need to be added with new arbitrary couplings. After 
this is done one continues to determine higher point amplitudes by recursion as usual. Because of 
this, one could have discovered the class ([^ by working solely with the scattering amplitudes, without 
recourse to any Lagrangian, as is the case with Yang-Mills theory. 
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Our other remark is that such constructibility is certainly not the property of a general effective 
field theory for the gauge field. Such a theory contains higher derivative interactions, and this prevents 
the amplitudes from having strong enough falloff properties under the so-called Britto-Cachazo-Feng- 
Witten (BCFW) shifts [6]. As the result, the BCFW recursion does not work for amplitudes in an 
effective field theory of the gauge field, at least not directly. The (generalised) constructibility of Q, 
in addition to their nice behaviour under renormalisation, is our other main motivation for their study. 

Yet another motivation for considering the class of theories Q is their relation to gravity. In [7] 
a certain infinite-parametric family of “deformations” of General Relativity was introduced by one 
of the present authors. In [5] they were recast as a family of diffeomorphism invariant SO (3) gauge 
theories, with GR being a member of this family [9]. In [inj some graviton scattering amplitudes were 
computed using this “pure connection” formalism. In this paper we shall see that deformations of 
General Relativity [8] continue to exhibit the 

Gravity = (GaugeTheory)^ (11) 

pattern known to relate GR and YM, with appropriate squares of the deformed YM scattering am¬ 
plitudes giving deformed gravity amplitudes. The fact that the double copy structure of gravity also 
extends to (some) higher order terms in the gauge theory and gravity Lagrangians is not new; it has 
been observed and studied in particular in m- Thus, our observations in this regard confirm the 
pattern noticed earlier. 

1.4 Unitarity and Interpretation 

Our final introductory remark will likely make some readers less interested in the whole story, but it 
is difficult to avoid mentioning from the outset. In Euclidean (or split) signature the self-dual part of 
the field strength is real. So, at the very least the theories we study make sense as Euclidean (split) 
signature field theories with intriguing properties. 

In Lorentzian signature the self-dual part of the field strength is complex. This makes the La¬ 
grangians Q complex. This immediately raises the question about unitarity and more general inter¬ 
pretation of the resulting quantum theory. In the next section we will give a quick Hamiltonian analysis 
of 0, and see that the corresponding Hamiltonians are not Hermitian, but instead parity-time (PT)- 
symmetric. Hence one is in the domain developed by the proponents of the HT-symmetric quantum 
mechanics m- As is explained in e.g. one way to give physical interpretation to such systems is 
to find a new positive definite inner product that makes the Hamiltonian a Hermitian operator. As 
the study of concrete PT-symmetric quantum system shows, see e.g. m, the equivalent Hermitian 
Hamiltonian can be expected to be extremely complicated, and so in practice for many purposes, e.g. 
determination of the energy levels, it is best to study the original non-Hermitian system. However, the 
non-trivial positive definite inner product that makes the Hamiltonian Hermitian plays an important 
role in the physical interpretation of the PT-symmetric quantum system. 

In this paper we make no attempt to determine (or establish existence of) the inner product that 
makes our quantum theory unitary. We shall study it as is, working with the obvious inner product in 
which the S-matrix is not unitary. The corresponding unitary S-matrix is then to be found from the 
non-unitary one by an appropriate similarity transformation. Alternatively, our scattering amplitude 
calculations can be interpreted as those for the split signature, where everything is real, but no physical 
interpretation is available. 

Thus, we make no claim that the quantities computed in this paper are ones that can be compared 
to some experimental measurements. Our interest in the class of theories Q is not because we want 
to use one of such Lagrangians to describe what happens in an experiment, at least not yet. Eor now 
we are interested in these theories because of the striking properties that they exhibit. Thus, our 
main goal here is to characterise the theories 0 in as much details as possible. Determination of 
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an appropriate physical interpretation of Q, if any, is only possible when their properties are well 
understood, and so is left to the future. 

In addition to non-unitarity discussed above, our theories are also plagued by more conventional 
unitarity problems arising in any non-renormalisable theory. This is the tree-level unitarity violation 
arising because in theories with negative mass dimension couplings the tree level scattering amplitudes 
at sufficiently high energies become larger than unity. This problem can only be cured (if at all) by an 
ultraviolet completion of the theory. The main objective of this paper can be paraphrased as to make 
a step in the direction of finding this UV completion for Q, by studying the arising renormalisation 
group flow and hoping for a UV completion in the sense of asymptotic safety of Weinberg [T] . 

1.5 Organisation of the paper 

We start in section with the Hamiltonian analysis of our theories. This illustrates what kind of 
complex Lagrangians we deal with. We also state the mode decomposition for the gauge field here, as 
well as fix our conventions, in particular for the helicity spinors. Section then describes two different 
gauge-fixings, one useful for loop calculations and the other useful for colour-ordered Feynman rules. 
We also derive the central “collapsing” property of our perturbation theory here, and evaluate the 
cubic vertices on-shell. In section [4] we discuss how BCFW recursion relations work for our theories. 
The main point here is that the amplitudes with at least one negative and one positive helicity gluons 
continue to be determinable by the usual BCWF recursion. The exception is the amplitudes with 
all plus helicity gluons. For these we need a more involved shift on all momenta. Section uses the 
BCFW recursions to compute the 4-gluon scattering amplitudes, and section discusses their double 
copy structure. In section we analyse the power counting for our theories, and illustrate it on an 
example of a divergent self-energy diagram. Section computes the triangle diagrams necessary to 
extract the beta-function of the new coupling. We conclude with a discussion. 


2 Hamiltonian formulation and the mode decomposition 

In this section we perform a quick Hamiltonian analysis of our theories. One justification for this 
exercise is to see what kind of complex Hamiltonians result. Another justification is to fix conventions 
and then state the mode decomposition of the gauge field. In the next section we use this mode 
decomposition to obtain the Feynman rules for computing the S-matrix. 


2.1 First order formulation 

To perform the Hamiltonian analysis of ([^ it is most useful to start with an equivalent but first order 
in derivatives formulation. Introducing an auxiliary field we write 

£ = - M*V{B/M^), (12) 

where V{B) is a gauge- and Lorentz- invariant function of the auxiliary field. The field equation for 
the auxiliary field is F^^ = M^(W)^^, where V is the matrix of first derivatives of V. It is clear 
that V{B/M‘^) is simply the Legendre transform of f{F/M‘^). 


2.2 Hamiltonian formulation 


We now perform the 3 -|- 1 decomposition of (12). In Lorentzian signature the relevant formulas are 

(13) 


T?a 

^MN 


= (2iFo“ - e*^'=FA)4 


jk) MN ’ 
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where i,j,k are the spatial indices and are multiples of Pauli matrices hxed by the following 

algebra they satisfy 


1 1 .. 

__ ^AC _j_ _ ^jkrpkAC 

2 ^/2 


(14) 


We also have 


FS, = A^-d,Al (15) 

where di is the covariant derivative with respect to the spatial connection. Introducing the new 
auxiliary field so that 

^MN = 

we can rewrite the Lagrangian as 

c = - n, (17) 

where 

= l^ijkFjk ( 18 ) 

is the magnetic field and the Hamiltonian is 

n = + M^V{P/M^). (19) 


We now see that the held is the canonically conjugate momentum to the spatial gauge heldH“. 
For Yang-Mills theory we have 

M^Vym{B/M^) = = ^(P“)2. (20) 

Thus, in this case we can rewrite the Hamiltonian as 

-Hym = y f7^“* - ^77“j +^(77“)2. (21) 

We recognise the canonical transformation of the usual Yang-Mills Hamiltonian. It can be checked 
that the above shift of the momentum variable is indeed a legitimate canonical transformation in that 
the symplectic form is unchanged. Thus, the correct reality conditions in this case are that the spatial 
connection, and thus the magnetic held are real, and the combination 


^ai ._ pdi _ _ pQ 


( 22 ) 


is real. This renders the Hamiltonian real, and theory unitary. 
For the theory Q deformed by a cubic term we get 


n = Sl ^ ^(77-)2 - 


a^^^2 _ ^abc^ijk 

3M2 


TT'CL I TJC 

' 2 




El + 


(23) 


plus higher order terms. Thus, if we require E and H to be real, the Hamiltonian is not Hermitian. 

(24) 


However, if a G M the Hamiltonian (23) does exhibit an anti-linear discrete symmetry 

pai ■ 


-E° 
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this followed by the complex conjugation. It is clear that this is the case for any V{B/M‘^), or 
equivalently f{F/M‘^) that contains only real coefficients when expanded into powers of its argument. 
In what follows we shall assume that the functions /, V satisfy this requirement. 


A remark about (23) is that this way of writing the effect of deformation makes it clear that at 


least the cubic deformation is only possible for the non-abelian Yang-Mills theory. Indeed, in the case 
of Maxwell we would not have an additional Lie-algebra index on the fields and the vector 


product present in (23) would simply be zero. So, for a U(l) gauge held the simplest new interaction 


is quartic in the held strength. 


2.3 Pseudo-Hermitian Hamiltonians 

As is explained in the literature on PT-symmetric quantum mechanics, see e.g. |13j . it is in the presence 
of an anti-linear symmetry that one may hope to make sense of the arising quantum mechanical system. 
In practice, this requires hnding a metric operator r] satisfying 


= rjHr] 


-1 


(25) 


where the dagger is dehned with respect to the ’’obvious” inner product, i.e. the one with respect 
to which the operator H is not Hermitian. The operators for which such r/ exists are called pseudo- 
Hermitian. One then dehnes the new inner product 




(26) 


with respect to which the Hamiltonian is Hermitian. If the new inner product happens to be positive- 
definite, it defines a new Hilbert space where the original Hamiltonian acts as a Hermitian operator. 
Mostafazadeh, see m and references therein, has shown that a necessary and sufficient condition for 
an existence of r] for a Hamiltonian H with a complete set of eigenvectors is that it commutes with 
an invertible anti-linear operator. 


In our case this operator is that given by (24) followed by the complex conjugation. Thus, one may 


hope for an existence of an appropriate metric operator r/. This would give one way of giving physical 
interpretation to our theories. However, as we have already explained in the Introduction, our main 
aim here is not to develop such a physical interpretation but instead study (formal for now) properties 
of theories 0. Thus, even though all calculations below are performed in Lorentzian signature and 
the arising S-matrix is not unitary in the obvious inner product, our treatment can be justified by 
going e.g. to the split signature. In this case there is still a scattering theory of massless particles 
that can be set up, but no unitarity to be concerned about. Of course our treatment of the scattering 
theory would be more compelling if we also found an inner product that makes everything unitary. 
But this is a difficult problem that likely requires a very different sort of considerations from the ones 
given below. We hope to return to this problem in the future. 


2.4 Mode decomposition 

We finish this section by stating the mode decomposition for the gauge field. We first note that all new 
terms in the Lagrangian only contribute to the interactions, and the kinetic term for the gauge-field is 
unchanged from what it is in Yang-Mills theory. Thus, the free Hamiltonian continues to be given by 

-Hfree = + ^(P“)^ (27) 

where we have only included into the Hamiltonian the part of the magnetic field linear in the gauge- 
field. The Hamilton-Jacobi equations can be written as 





where the operator ed acting on any Lie-algebra valued spatial vector is defined as 

{edXf^ := 

It is an easy exercise to check 

(e5)2x“* = d^djX'^^ - d^djX'^\ 


( 29 ) 


(30) 


Thus, on transverse vectors it reduces to minus the Laplacian. Therefore, as expected, the two 
equations (28) together imply that the transverse = 0 part of the connection satishes the wave 
equation. 

We can now write the mode decomposition for the transverse part of Af 

d?k 


A<l = g 


(27r)32a;fc 


+ £+(o+)“e'‘"’ + 


—ikx 


(31) 


The factor of g in front will of course later be absorbed into the gauge field perturbation. Our signature 
is —, +, -)-, + so that kx = k'^Xi — ut. For our massless gluons LOk = \k\. The vectors are the two 
circular polarisation helicities. They are chosen to satisfy 


I-- 






ICi 




(32) 


and 


4 = 4-)' 


(33) 


The condition (32) guarantees that only the positive polarisation gives a non-zero momentum vector 

pai 

The decomposition (^3l| together with the equal time commutation relations [A^ (x), (y)] = 


{x — y), where Iljj is the transverse projector, implies the correct commutation relations 

[i4r, )^'] = {2nf2uJkS'^»6^{k - p) (34) 

provided the following normalisation condition is satished 


e+e7 + e-e+ = H 


* j 


* j 


Hj- 


(35) 


Thus, in a coordinate system in the momentum space in which the z-direction is taken to be along 
the direction of the momentum vector ki, the polarisation vectors are given by 


4 = 


(36) 


In particular = I. 


2.5 Mode decomposition in the spinor form 

The last equation we need to set up the theory is the mode decomposition in the spinor form. Writing 
^MN ~ where the matrices have already been introduced above, and then following 

the rules for converting into a space-time form described in subsection 6.10 of [TO], we get 


A 


MM' — 9 


d^k 


(27r)32a;fc 






-ikx I - 

^ ^MM' 


( 4 )'“ 


—ikx 


(37) 
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with 


^MM> 


= 1 - 




{qk) 


-MM' 


= 1 - 


. kMQM' 


[q k] 


(38) 


Here (A^) := A^^/i^JA//] := are the two spinor contractions, and qM,qM' are the reference 

spinors whose presence reflects the gauge transformation freedom. It is easy to see that the required 


normalisation condition is satisfied ~ 


We now have the usual rules of the Lehmann-Symanzik-Zimmermann (LSZ) reduction, with the 


understanding that the helicity spinors (38) must be inserted into the amputated Feynman diagrams 
to extract the scattering amplitudes. 


3 Gauge-fixing and the Feynman rules 

There are two useful gauge-fixings that we will consider that lead to two different versions of the 
Feynman rules. One of them is more useful for one-loop calculations, the other is useful to derive 
the colour ordered Feynman rules. Since the new terms that we added to the Lagrangian only affect 
interactions, the gauge-fixing(s) that we need to introduce are the same as for YM theory. The only 
slight novelty here is that we work with a chiral formulation of YM, so even the resulting YM Feynman 
rules are slightly different from the usual ones. In particular, spinor notations will be crucial for writing 
all formulas. 


3.1 Feynman gange 

Absorbing the YM coupling into the gauge field perturbation, the kinetic term for the Lagrangian ([^ 
reads 

. (39) 

Using the magic of spinor identities we can rewrite the second term in brackets as 

Onm'Am^ = dMM'A%^ -I- eMNd^ m'A%^ . (40) 


This formula is easily checked by multiplying both sides with 
raising and lowering of the spinor indices = A^, X^emn - 

indices. Thus, we can rewrite (39) as 


We also use the following rules for 
Xjsf, and the same for the primed 




(41) 


where our convention for the metric contraction is m'Vm^' ■ In the first term here there is 

a contraction of the partial derivatives in an unprimed spinor index. Given that the partial derivatives 
commute we have 

N' = (42) 

This formula is easily checked by multiplying both sides with , Thus, overall, modulo surface 
terms 

= ( 43 ) 

as expected. 
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We now add a gauge-fixing term that removes the last term, and add the ghost Lagrangian 

^ghost = (44) 

where is the covariant derivative with respect to the gauge field, and c“,c“ are the ghost and 
antighost fields. Overall the gauge-fixing term and the ghost Lagrangian are exactly as in the usual 
YM theory. 


3.2 The Feynman gange Feynman rules 


The momentum space propagator for the gauge field follows directly from the gauge-fixed kinetic term 
and reads 


Hi uj 


1 

ik^ 




(45) 


The extra minus sign here is due to our convention for the metric Hmm'nn' = —(-mn^m'N'- 

To derive the vertices we first note that the derivative of the gauge field only appears in the 
combination (10). For reasons that will become clear below, it is then convenient to represent the 
cubic YM vertex as a sum of three terms, where in each of the terms the derivative is applied to a 
different external leg. We denote the leg where the derivative acts with a dot. The resulting Feynman 
rule is then 



We have defined 


{kA{k)f^^ := 2k^^A^^'^^'{k) 


(46) 


(47) 


(48) 


and assumed that all the momenta are incoming, and A{k), A(p), A{q) are the placeholders. These 
can later be either replaced with external polarisations, or Wick contracted with another placeholder 
to form the propagator. 

The YM quartic vertex reads 


k a 


p e 



I b 


q d 


= + 2 terms (49) 


The additional terms that we did not write are analogous, just with different pairs of placeholder 
gauge fields contracted into the structure constants. 
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The cubic and quartic vertices following from the new cubic term (second term in Q) are as 
follows. First, there is a completely symmetric cubic interaction 



= ^J'^^%kA{k))\^{pA{p)%^{qA{qW 


Then there is the following quartic vertex 

k a ^ ^ lb 


p e 


X 


+ 5 terms 


q d 


where 


k a 


p e 


X 


I b 


q d 


^9 rabc fcde / 




(50) 

(51) 


(52) 


(53) 


The new term in ([^ produces also quintic and sextic interaction vertices, but we will not need them 
here. 


3.3 Central identity 

The benefit of using the dot notation consists in the fact that when two dots appear on the same 
internal line, this line collapses to a point. Indeed, consider the following Wick contraction 


^ {kA°')AB {kA^)cD ) = ( ) 

= —4:5°‘^k(^AA'^^^'"£B){Dkc)C' = -rb°'^^A(C^\B\D)- ( 54 ) 


We thus see that the k"^ in the denominator gets cancelled by arising in the numerator. The way this 
happens is instructive. To understand this we need to trace the contraction of the primed indices only. 
These are contracted in each of the two factors of (kA). Then the contraction of the two placeholder 
^d’s makes the primed indices of the two factors of k contract. This produces a factor of times an 
e of the unprimed indices. It is this k'^ appearing in the numerator that cancels the propagator. 

The identity (54) can be represented graphically as follows 



(55) 


where an additional minus sign comes from the fact that on the internal line one of the momenta is 
incoming and the other is outgoing. The dots represent any lines that may or may not be present on 
the left and the right. The lines that are external in this diagram may also contain derivatives, i.e. 
other dots. Note that the effective vertex arising needs to be worked out on a case-by-case basis, see 
examples of this below. 
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The identity (54) clearly holds due to the very special way (10) that the derivatives of the gauge-field 
appear in the interaction vertices. Clearly, no such identity would hold if just the anti-symmetrised 
derivative would be involved, as in the usual Yang-Mills Feynman rules. This is the benefit of 

the “chiral” formulation considered here. As we shall soon see, it is the identity (55) that is at the 
heart of the improved properties of our theories as compared to a general effective field theory. 

The identity (54) is not new; it has been noticed in particular in jlSj, see (5.10) of this paper and 
in m see (4.15). The “dot” notation introduced here, however, appears to be new. This notation 
is particularly useful for analysis of the arising diagrams, as well as their cancellations. Thus, using 
our notation it is for example very easy to verify that the only purpose of existence of the YM theory 
quartic vertex is to cancel the YM cubic vertex diagrams where the two dots happen to appear on 
the same internal line. Because of this fact, the usual YM theory can be recast into a theory with 
only cubic vertices (46) as well as the rule that there is no more than one dot on any line. This fact 
appears to be important, and will be explored elsewhere. 


3.4 Gervais-Neveu gauge 

To derive the colour-ordered Feynman rules convenient for the scattering amplitude calculations, we 
introduce a different gauge. We basically follow m Chapter 79, with the only difference being that 
spinor indices are used instead of space-time, and the original YM Lagrangian is chiral. 

Let us introduce the matrix-valued generators r“ satisfying 

Tr(r“r'') = [r“,r'’] = i\/2/“^'=TL (56) 

We define the matrix-valued field 


Amm' 

The self-dual part of the matrix-valued curvature is then 

9 


FmN — 

and the chiral Yang-Mills Lagrangian is 


+ +A 


NM'^M 


A 

Aif - 


(57) 


(58) 


Fym = -^Tr(F^^FMN)- 


(59) 


The gauge-fixing term is supposed to make the kinetic term for the gauge field non-degenerate. 
Previously this has been achieved by adding —(1/2)(9^A“)^ to the Lagrangian. This only modifies 
the kinetic term, leaving the interaction vertices unchanged. However, it is possible to add terms 
non-linear in the gauge field to the gauge-fixing term so as to modify the interactions as well. A 
convenient choice is 




H = d>^A^ + 

iv 2 


(60) 


Indeed, by squaring (58) both the cubic and the quartic vertices contain two terms, coming from the 
symmetrisation of the MN indices. The choice (60) cancels one of the two terms in both cubic and 
the quartic vertex. Let us see this. Using the same trick as in (40) we can rewrite the self-dual part 
of the curvature as 


Fmn = 2 ( Omm'An^ + -tj^Amm'An^ 


+ ^mnH. 


(61) 
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Then 


-^^YM = —Tf ( dMM'^N^ + 




and hence 


Cym + C^.L = -\T,{d^^A'^d^A,) + 6^1 + 4 m 


(62) 


(63) 


with 


■^YM = \gV2Av{d’^A mn'^n^ ), 


C 


( 4 ) 


,ym = Amn'An^^ ). (64) 


We will also need the new vertex from Q in the same matrix notation. Writing in terms of gauge 
fields defined with one factor of coupling constant absorbed, we have 


Tnftw — i 


ag 


-Tt{Fa^Fb^Fc^). 


3V2M' 


(65) 


3.5 Possibilities at higher order 

This representation of the new interaction also suggests a way to encode and classify the possibilities 


at higher powers of the curvature. Indeed, if {T°‘)i is an N x N matrix, it is clear that in (65) the 
curvature Fa^ is thought of as a 2N x 2N matrix 


Fa^Braj _ 


( 66 ) 


SO that the interaction (65) is just a multiple of the trace of the cube of F. Thus, at the order quartic 


in curvature we can have a single trace interaction Tr(T^), as well as the double trace interaction 
Tr(F^)Tr(.F^). This way of classifying the new interactions can be clearly extended to higher orders. 


3.6 Colour-ordered Feynman rules 

As we just saw, at quartic order and higher there are also multi-trace interactions in our class of 
theories. Because of this, even at tree level the scattering amplitudes for our theories no longer follow 
the simple single-trace pattern. This is not a cause of concern for us here, because this more involved 
structure will not show up in the low point amplitudes we compute. But one should keep in mind 
that the colour ordered rules for our theories are more involved than in YM because of the presence 
of multi-trace interactions. 

At the cubic interaction level all interactions are single-trace, and we can easily deduce the colour- 
ordered Feynman rules from (64) and (|65l) . We will be determining the tree-level scattering amplitudes 


using the recursion relations, and so we will only need the Feynman rules for these cubic vertices. We 
the colour stripped off, the YM cubic vertex is again represented as a sum of three terms, with the 
derivative in each term acting on a different leg. We again use placeholders, but these are now colour 
stripped. We have for the YM cubic vertex 


- igV2kM'A 


NM' 


{k)AMN'ip)AN^ {q) + 2 terms 


and for the new vertex 


ag 

‘\/2M2' 


- 1 ^ {kA{k))A^{pA{p))B^{qA{q))c" 


(67) 


( 68 ) 


where 2k^ m'A^^ (k) = {kA{k))^^. As before, the propagator is a factor of l/i/c^ times a set of 
Kronecker delta’s as well as spinor e’s. All momenta are assumed incoming. 
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Note that the colour-ordered Feynman rules (67) and (68) also admit a representation analogous 


to (46) and (50). In other words, one can still think in terms of the derivatives acting on legs of a 


diagram in terms of dots. One still has the key property (55) that two dots on the same internal line 


collapse the line. This will be important in what follows to estimate the high energy behaviour of the 
diagrams. 


3.7 Colour-ordered cubic vertices on shell 

On-shell cubic vertices in any theory of massless particles vanish because of the special kinematics 
for three real null momenta satisfying the momentum conservation. However, as has become the 
common practice in the literature on scattering amplitudes, it is convenient to analytically continue 
the momenta to complex values in a certain way. One then obtains non-zero answers for the on-shell 
3-point amplitudes. These then play a crucial role in the proof of constructibility, as they serve as 
seeds for the recursion that is used to construct all higher point amplitudes. 


Thus, we now evaluate the vertices (67) and (68) on shell (with momenta continued to complex 
values), as this is the input of the recursion relation calculations. We first note that 


{ke 


-\MN 


= 0 , 


{ke+) 


MN 


= -2ifc“A: 


MuN 


(69) 


(70) 


Hence the new vertex is only non-zero for the all plus configuration 

Att+ = ^^^(12)(2 3)(31). 

Yang-Mills vertex, on the other hand, is zero for both the all plus and all-minus configurations. For the 
all minus configuration this is clear, as there must be at least one plus helicity so that the derivative 
term does not give zero. For the all plus configuration, there is a contraction of the primed indices of 
the two gauge field placeholders that are without the derivative, and this is what makes the amplitude 
zero. 

Let us compute YM vertex for the single plus two minus configuration. This is easy, because the 
plus helicity state should be inserted where the derivative acts. We get 

-^YM - 9^q) (3 q) ^ ^ [i 2] [31]’ 

where we have used the momentum conservation in the last step. 

The opposite helicity configuration is a bit more complicated because there are two terms. We get 


A — 

•^YM 




We now use the momentum conservation to get 


. 4 -++ - qV2^^ 


(73) 


which is just the complex conjugate of © as it should be. 


4 Recursion relations 

In this section we explain how the tree level scattering amplitudes can be computed using recursion 
relations. First, the maximally helicity violating (MHV) amplitudes, now defined as those with just 
two plus helicities and the rest minus are unchanged from the YM case. Then the amplitudes with 
at least one minus can be computed using the usual BCFW recursion where a pair of momenta are 
shifted. The all plus amplitudes, which are now non-zero, are computed using the Risager’s shift 
BCFW recursion. 
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4.1 Amplitudes with all negative helicity gluons 


It is easy to see that amplitudes for all negative helicity gluons, or just a single positive helicity 
configuration, continue to be zero for our theories as they are in the case of YM. Indeed, as usual, we 
choose all reference spinors of negative helicity gluons to be the same. In this case, the only hope to 
contract these helicity spinors and get a non-zero result is to use the momenta in the numerators of 
the diagrams, coming from the derivatives present in the vertices. This is best achieved by having as 
many derivatives as possible in a diagram. And the maximum number of derivatives in a diagram is 
achieved by the 3-valent graphs. 


It first looks like the new vertex (68) should change the situation drastically, because it contains 


three derivatives instead of single derivative in the YM vertex case. But because of the identity (55) 


there cannot be too many copies of this vertex in a diagram. Indeed, the identity (55) implies that 


there can be at most a single dot on every internal line of a diagram. Consider now a 3-valent diagram 
for n negative helicity gravitons. There are n — 3 internal lines. Note that the dot cannot appear 
on the external lines, because it will then give zero applied to a negative helicity graviton. So, the 
maximal total number of derivatives in a diagram for all negative helicity gluons is n — 3, which is not 
enough to absorb all n reference spinors. This shows that all negative helicity amplitudes are zero. 

A similar argument is at work for the configuration with at most single positive helicity gluon. 
We can then choose the reference spinors of negative helicity gluons to be the unprimed momentum 
spinor of the positive helicity gluon. Again, we have n unprimed spinors to contract, while avoiding 
contracting them to each other. Now that we have one positive helicity gluon, we can have one 
derivative on an external leg, raising the total number of possible derivatives in a 3-valent diagram to 
n — 2. This is still not enough to contract all n reference spinors, hence giving a zero result for such 
amplitudes. Below we will also treat this case using the BCFW recursion. 


4.2 Amplitudes with at least one negative and one positive helicity gluon 

Let us now consider more general amplitudes. Let us assume that there is at least one negative 

helicity gluon, and at least one positive. The case not covered is thus just that of all positive helicity 

configuration. This will be treated separately. 

We would now like to analyse the behaviour of amplitudes under the BCFW shift, where the shift 
is applied to a negative-positive pair of gluons. Let us label the gluons so that a chosen negative 
helicity gluon has label 1, and the chosen positive gluon label n. We then perform the BCFW shift [6] 

I —i = I -|- z n, n ^ n = n — zl'. (74) 


Under this shift the helicity spinors (38) of both gluons change and behave as 1/z for large complex 
parameters 2 ;. 

We can now count the powers of 2 ; appearing in a generic Feynman diagram shifted as in (74). We 
are interested in the large 2 ; behaviour, and would like to show that the amplitudes go as I/ 2 ; in this 
limit. Thus, as in the previous subsection, we need to be concerned about the derivatives present in 
the vertices. The worst case scenario is when there are only cubic vertices present. However, as we 
already used in the previous section, in this case there can at most be n — 3 derivatives on the internal 
lines. This equals to the number of propagators. Each propagator on a curve connecting vertices 1 
and n will contribute a factor of I/ 2 ;, and the derivatives on these internal lines will each contribute 
a factor of 2 ;. Thus, there will be a cancellation of factors of 2 ; from the internal lines. The last thing 
to note is that both external lines connecting to 1 and n cannot contain a derivative. Indeed, as one 
of the helicities is negative there cannot be a derivative on that external line. Thus, we see that the 
helicity spinors each give a 1/z factor, while there is at most a single derivative on the external line, 
giving another factor of 2 ;. The above discussion can be summarised in a formula as follows. Assuming 
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1 - 


n+ 



Figure 1: Standard BCFW shifting on the pair (1 ,n~^) 


the number of internal lines connecting the shifted pair (1 , n"*") is I, we get the following z-behaviour: 


- z X z X \ - ] = - 


(75) 


Overall, we see the desired 1/z behaviour of such amplitudes. An illustration of this shifting is provided 
by Fig 1. 

It is also not hard to deduce the desired 1/z behaviour of the amplitudes from an argument 
analogous to the one used in [15]. This is the background field method argument, where the action is 
expanded to second order in the “hard” field, on a background of “soft” particles. To obtain the 1/z 
behaviour one should just note that the interaction of the hard perturbation with the soft background 
contains at most a single derivative of the hard field. This is trivial to see in the YM case. In our 
“deformed” case, there is also a second-derivative interaction given by 

a pa ipb i9^)MNidA)pQ, (76) 


where the derivative of the function / is evaluated at the background. The other, subleading interaction 
terms contains at most a single derivative. It seems that we now have a dangerous second-derivative 
interaction of the hard and soft fields. What saves the game is that for the helicity configuration —h 
one of the hard fields is necessarily a negative helicity field, for which the object {dA)^j^ vanishes. 
So, for the —h hard field configuration, the interaction with the soft field starts with a first derivative 
term, which is what is needed to get the desired 1/z behaviour. 

The conclusion is then that the amplitudes with at least one minus and one plus gluon behave 
as, 1/z under the BCFW shift (74). Because of this they can be determined by the BCFW recursion 
relation, with only on-shell vertices needed to compute the amplitudes. However, to run the BCFW 
recursion, generically we also need to know the all plus amplitudes. 


4.3 MHV amplitudes 

We define the MHV amplitudes as those with exactly two positive helicity gluons, and an arbitrary 
number of the negative helicity particles. These are the simplest non-vanishing amplitudes. To 
determine these amplitudes we do not need to know the all plus amplitudes, and so we can discuss 
them already now. 

Let us first note that we can use the BCFW recursion to show that amplitudes with just one plus 
helicity gluon continue to be zero. Indeed, such amplitudes are cut by BCFW into two copies of lower 
point amplitudes of the same kind 

M(l-, 2-,..., n+) = 2-,..., fe-, P+)^A{-P-, /c + 1",..., n+) = 0, 

k 
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and thus are zero by induction. 

It is now easy to see that the MHV amplitudes continue to be given by the same expressions as 
arise in Yang-Mills theory. This follows from the applicability of the BCFW recursion relations in 
this case, as there is clearly at least one minus and at least one plus helicity gluon. Given that the 
all minus and all-except-one minus amplitudes are all zero, the on-shell relation for MHV amplitudes 
prescribes 

... ,i+,... ,n - l",n+) = A{i~,2~, ... ,i+,..., re - l“,re+) 

+^(1-, 2-, P+)^^(-P-,3-,... ,i+,... ,re+). 

The first line here vanishes at the pole, while the second line gives the same contribution as in the 
original Yang-Mills theory. So, one is never using any new amplitudes in this recursion, and the arising 
MHV amplitudes are just those of Yang-Mills theory. 


4.4 All plus amplitudes 


Our above arguments based on the BCFW shift (74) are not applicable to the all plus helicity con¬ 
figurations. And indeed, we should not expect these amplitudes to decay under the BCFW shifts. 
Consider for example the (74), but now applied to two positive helicity gluons. It is clear that under 
this shift the basic 3-point amplitude (70) behaves as 0{z) for large z. Thus, the simple BCFW shift 
acting on just a pair of gluons is certainly not going to work for the all plus amplitudes. 

However, there is another possibility, first used in work by Risager m. where all gluons are shifted. 
Thus, let us shift the primed momentum spinors of all the positive helicity particles as 


i' = - 1 - zfiq', ^ ifi = 0. (77) 

i 

In other words, we shift each primed momentum spinor in the direction of the reference spinor used 
in the helicity spinors. This is done in such a way that the momentum conservation is respected. An 
explicit form for the coefficients /* will be given below at four points. 

We will see that the shift works for our purposes, because it can be used to determine the parts 
of the re-all plus amplitude that arise from lower order amplitudes. However, generically at re-points 
there will also be parts of the amplitude that cannot be determined by the recursion, the so-called 
BCFW boundary terms. However the sources of boundary terms can be readily identified. Indeed, 
assuming for simplicity that we deal with the theory with just a single new interaction (§, it is easy 
to see that the only diagrams not exhibiting the l/z behaviour are those constructed solely from new 
cubic vertices. The contribution to BCFW boundary is then simply given by the sum of diagrams 
that completely collapse 





This sum can be evaluated explicitly, as the terms that appear are just spinor contractions {ij). An 
example at four points will be considered below. 

Alternatively, one can note that the boundary terms that appear are just multiples of the higher- 
valent vertices that anyway have to be included into the theory. So, one learns that the parts of 


18 





the n-point all plus amplitude that are not determined by the BCFW are just the n-valent vertices 
evaluated on shell. These parts come with new coupling constants and should simply be added to 
the those parts that are determined by the BCFW recursion with Risager’s shift to obtain the full 
amplitude. 

To see that the rest of the diagrams indeed converge under the all leg shift (77), let us apply 
the collapsing rule until there are no internal lines to collapse. If no internal lines are left, we are in 
the situation of the boundary term just considered. So, we are left with diagrams with “effective” 
higher-valent vertices, and internal lines with just a single dot or no dots. The lines with a single dot 
are dangerous, because there is a factor of 1/z coming from the denominator, and a possible factor of z 
from the derivative. Lines without dots are not dangerous, as they immediately give a 1/^; behaviour 
to the diagram. Thus, we can concentrate on dangerous diagrams where all internal lines carry a 
single dot. 

To see the 2 ;-behaviour of these diagrams, let us consider a line with a single dot on it. This dot 
is on one end of the line, and so there is no dot on the other end. Thus, the end of the line that does 
not have a dot is contracted into some effective vertex of the type 



Note that there is necessarily two (or more generally even number of) legs without dots, as otherwise 
the primed index of the placeholder connections will have nowhere to contract. We now follow the 
other leg of this vertex where there is no dot. If there is more than a single pair of undotted legs, we 
follow the leg to which the primed index of the placeholder for the line we came from is contracted to. 
We concluded that the dangerous diagrams are those where every internal line has one dot, so there 
is a dot on the other end of this second line. 

What we now have is a pair of internal lines, each with a single dot on them, connecting at an 
effective vertex. What is important for us is that in this effective vertex the primed indices of the 
two gauge field placeholders contract. This means that the numerator arising from the dot on the 
first line is contracted to the numerator on the second line in their primed indices. As the result, the 
would be dangerous behaviour of the product does not arise, because it is necessarily proportional 
to [q q] = 0. This establishes that there are no other boundary terms, and everything apart from the 
fully collapsed diagrams falls off as l/z or even faster as z —)■ oo. So, we can use Risager’s shift to 
determine the “complicated” parts of the n-point all plus amplitudes. 

5 4-point examples 

The purpose of this section is to illustrate the general statements of the previous section on the 
examples of 4-point amplitudes. Some amplitudes with Fgj^ interaction have been previously computed 
in the literature; see in particular m- So, there is no claim for novelty in this section, and it is given 
more for illustration on the types of amplitudes that arise. 

5.1 (-1-+) MHV amplitude 

As we have already discussed, there are no all minus and-h amplitudes, so the first one to 

consider is the MHV amplitude (-!■+)■ We expect that it is given by the same expression as in 

the usual YM theory. 
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From the Feynman diagram calculation perspective it may seem that the following two new dia¬ 
grams may contribute to this amplitude 



However, it can be checked that they precisely cancel each other in view of the collapsing rule. So, 
the MHV amplitude is completely determined by the original Yang-Mills diagrams. 

The same conclusion can also be derived from the on-shell relation perspective. As we discussed 
in the previous section, we can apply the —h shift (74). The only possible contribution is then 


M(l-, 2-, 3+, 4+) = M(i-, 2-, P+)^M(-P-, 3+, 4+), 


(78) 


with = 2(1 2)[12]. Thus, the on-shell recursion relation of the MHV amplitude does not contain 
the new 3-point all plus amplitude and is therefore unchanged in the deformed theory. 


M(l-,2-,3+,4+) = {-i)g^ 


[ 12 ]^ 


[12] [23] [34] [41] 
and similarly for the non-adjacent MHV amplitude. 


= 


( 34 )^ 


(12)(23)(34)(41)’ 


(79) 


5.2 (— h -I--I-) amplitude 

The first 4-point amplitude that is affected by the deformation is M(l“, 2'’', 3^, d"*"). At first sight we 
have one diagram that may contribute to the BCFW boundary 


1 - 


2-h 



4 + 


3-h 


However as in the previous discussion this diagram behaves as 1/z because of the dots collapsing rule. 
Therefore, the amplitude is given by an on-shell relation as 

M(l-, 2+, 3+, 4+) = M(i-, 2+, P+)^M(-P-, 3+, 4+) + ^(1", 2+, p-)^M(-P+, 3+, 4+) (80) 

The first term is not contributing because M(l“, 2'’“, P"*") vanishes at the pole, whereas the second 
term gives 

.4(i-.2+.p-)j2^(-F+,3+,4+) = V29jF:|l_^2=V2^(P3)(34)(4P> (81) 


Note that this is different from the usual Yang-Mills theory where we do not have a non-vanishing all 


plus 3-point sub-amplitude. Substituting the on-shell value for internal momentum P = A:2((4 1)/(4 2)A:(-|- 
^ 2 ) we get 


M(l-,2+,3+,4+) 


Aag^ (23)(34) 
* M2 (12)(14) 



(82) 
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5.3 The all plus amplitude 


When all four external legs carry plus helicities, the all plus amplitude diverges under the standard 
BCFW shift (74). Let us instead consider the following generalisation of the Risager shift on all four 
external legs 


ki = kik'i +0 

k2 = k2k2 +zk2{M)ri' 
h = /csfcg +zA:3(42)r/' 
^4 = ^ 4^4 +zk4(^23')r]' 


+2/01(34)7/' +2/01(24)77' 
+0 + 2 / 02 ( 41 ) 77 ' 

+ 2 ^ 3 ( 41 ) 77 ' +0 

+ 2 ^ 4 ( 13 ) 77 ' + 2 / 04 ( 12 ) 77 ' 


+2 /oi(23)77' 
+ 2 / 02 ( 31 ) 77 ' 
2 ^ 3 ( 12 ) 77 ' 
+0 


(83) 


Let us also choose all reference spinors to be identical to the same 77 ' introduced to define the shiftings. 
Since only the anti-holomorphic (primed) spinor are shifted, the polarisations are unchanged. 

Then the diagrams consisting of two new cubic vertices is finite for large 2 and therefore contributes 
as the BCFW boundary term(s), see below. All other diagrams go to zero for large 2 , and can be 
calculated separately. Indeed, the s-channel residue is given by the following two MHV/(+ + +) 
combinations, which both converge under Risager shift (but diverge in the standard BCFW setting) 


A(i+,2+,p-)^A(-P+,3+,4+) 

A(i+,2+,P+)^A(-.p-,3+,4+) 




For brevity we denote the shiftings (83) taken on legs 1 and 2 as 


ki = k'l + zcirj, 
^2 = k '2 + z € 277 ' 


(84) 


(85) 


The on-shell condition at the s-channel pole allows us to write 

1 


P = 


Cl [ 277 ] - C 2 [irj] 


{ki [I 77 ] + /c 2 [ 277 ]) (ci/c 2 - C 2 A:'i) 


( 86 ) 


The s-channel residue can then be readily obtained by substituting P = ki [I 77 ] + k 2 [ 277 ] = —k^ [ 377 ] — 
/c 4 [ 477 ]. We get 


(-i) 




( 12 ) 


(4 + 4), 


[3 4] [1 77 ] [2 77 ] [3 77 ] [4 77 ] 
where we introduced the following anti-symmetric quantities: 

dj := [ni]{ij)[jv]. 

Repeating the same calculation for the t-channel, we find the total amplitude 

1 


(87) 


( 88 ) 


7l.mg«Zar(l + ,2+,3+,4+) = (- 7 )^"^ 


M2 [ 1 77 ][2 77][3 77] [4 ))J (m ("'2 + 4 ) + (4 + ci2) 


Using the momentum conservation this becomes 

4ac/2 


7l.*n3«iar(l + ,2+,3+,4+) = (-i) 


( 12 ) 


m2 [3 4] [1 77 ] [2 77 ] [3 77 ] [4 77 ] 


( 4 + 4 - 4 - 4 ) 


(89) 


21 













The momentum conservation implies that we have the following identities satisfied by the c’s 


E 


— 0 * 


( 90 ) 


Because of this we can rewrite the quantity in brackets in (89) as 


(Ci2 — Ci 4 )(ci 2 + C14) + (C34 — C 23 )(C 34 + C23) = (ci2 — Ci 4 )(—C13) + ( —C3i)(c34 + C23) = 
C 13 (C 34 + C23 - C12 + C14) = 2C13C42 = 2[lr/] [2r]\ [3??] [4r/] (l3)(42). 


(91) 


Thus, the singular part of the 4-point all plus amplitude is finally given by 

4l.*n5nZar(l + ,2+,3+,4+) = (-*) ^ ^(13) (42) (12). (92) 

The BCFW boundary term, on the other hand, can be derived from a direct calculation of the s 
and t-channel diagrams consisting of new cubic vertices only. In viewing of the dot collapsing rule, we 
can denote the sum of these two diagrams as 





The boundary term is given by contracting these diagrams with the positive helicity polarisations. 
For the first one we get 

i-^(12>{34)({13>{24) + {23){14)), 

and the second one gives 


\QoP‘ 

i-^{14){3 2)({13){42) + {43){12)). 


The sum is 


24bo?in(iar?y(l~*~) 2^, 3~*~, 4~*~) — i 


M4 


(13)^(24)^ - 2(12)(23)(34)(41) 


(93) 


It is clear that the boundary term in the all plus amplitude is of the form that can be obtained by 
evaluating the 4-vertex of type on shell. As we discussed previously, at the quartic level there are 


two possible interactions, one single-trace, one double-trace. The last term in (93) is precisely of the 


type that would be obtained by evaluating the single-trace quartic interaction on-shell. The first term 
is of the double-trace type, but here one must be careful to include all colour orderings to interpret 
the arising colour structure correctly. In any case, both contributions to the boundary term are of the 
type arising by evaluating the quartic vertex on-shell. 

One now has two options. One option is simply to interpret this boundary contribution as a shift 
of of the new couplings that get added to the theory at the 4-point level. Then one simply forgets 
about this boundary contribution, and adds a set of all plus amplitudes with arbitrary couplings. 

Another option, possibly more interesting, could be to require that the amplitudes are completely 
determinable by the recursion, with no new couplings to be added apart from the already existing a 
coupling. Thus, one can require that the new interactions of type that are added to the Lagrangian 
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are designed so that there is only a contribution to the 4-plus amplitude behaving as 1/z under the 
all momentum shift. In other words, one can decide to cancel the above boundary term with what 
comes from the interactions of type F^. If one continues to add new interactions at every order so as 
to have the I /2 falloff of the all plus amplitudes, one obtains a set of amplitudes described just by 
two couplings, the YM coupling g and the new coupling a. One can also decide to add a hnite set of 
higher-order couplings, and then request the Ijz behaviour of all plus amplitudes starting from some 
order n. This way one will obtain amplitudes parameterised by a finite number of couplings. All these 
options are possible due to the constructible character of amplitudes for our theories. 


6 The double copy structure 

The purpose of this section is to point out that the 4-gluon amplitudes obtained in the previous section 
continue to be related to the 4-graviton amplitudes for deformed GR known from [TO] via the double 
copy pattern. The fact that non-renormalisable gauge theory interactions continue to square to 
gravity amplitudes is not new, and it has been discussed in e.g. [20] and m- The latter reference 
also contains a discussion of the colour-kinematics duality, which we do not go into. What is new 
in the present paper is the purely chiral context where both the gauge theory and gravity sides are 
augmented by only the chiral half of the amplitudes. 


6.1 Graviton scattering amplitudes for deformations of GR 

In this paper we will only look at the relation between 4-particle amplitudes, leaving the general case 
to future study. 

The 4-graviton scattering amplitudes in a family of “deformations” of General Relativity studied 
in m are as follows. The all minus and one plus amplitudes continue to be zero as in GR. The 
-1-+ MHV amplitude is unchanged from what it is in GR, and is given by 
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(94) 


with the Planck mass being Mp = l/lGvrG. The —h -|-+ amplitude is non-zero in deformed GR, and 
is given by 


4M6L J (12)(13)(14) 


(95) 


where [3 is the parameter for a particular one-parameter family of deformations described in the 
Appendix in [10] . Finally, the all plus amplitude consists of two contributions 


. 9 / 3 ^ 


= i^(12>(13>(14)(23)(2 4>(34>|| 

(^(12)^(23)^(34)^(41)^ -h (12)^(24)^(43)^(31)^ -h (1 3)^(32)^(24)^(4 1)^) . 


(96) 


The first of these contributions comes from a Feynman diagram calculation and has s, t, u plane 
singularities. The second contribution is just obtained by evaluating the quartic vertex on-shell. The 
coupling constant appearing in front of the quadratic vertex does not have to be related to the cubic 
vertex couplings. Thus, the coefficient in front of the term in the second line is in general arbitrary. 
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6.2 The double copy structure of the MHV amplitude 

Here we review the usual gravity as square of YM story. A good reference for this is [2T] . At the level 
of the 4-point amplitudes we have 

M^Af—++(1, 2, 3,4) = i4^—++(1, 2, 3, 4)si2\a—++{1, 2,4, 3), (97) 

C gZ gZ 

where here si 2 = (1 2)[1 2], Given the above explicit form of the amplitudes, it is a simple exercise to 
verify this relation. We have a different sign in this formula as compared to one in EH, which is due 
to our different convention for the helicity spinors. 


6.3 The double copy structure of the new amplitudes 


Let us apply the double copy relation (97) to the —h ++ amplitudes. In this case the left hand-side 
gives 

/4a (3 2)(3 4) A ;io'rio|G° (42)(43> A _ 16a^ (23)=‘(24>=‘(34>^ 

VVM12)(14)' V VV2(12)(13)' V A/‘ ' ' {12)(13)(14) ’ 

apart from an extra minus sign that is likely due to our different helicity conventions, which starts to 
play a role when there is a different number of positive and negative helicity particles, and this equals 


to the left-hand-side in (97) with the gravitational amplitude given by (95) if we identify 


16 a^ 

M4 
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4A^' 


(98) 


Let us also work out the square of the (singular part of the) gauge theory -|- -|- amplitude. The 
left-hand-side of (1^ in this case gives 
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8 a (13)(24)(34)\ 
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( 12 ) [ 12 ] 


8 a (14)(2 3)(4 3)\ Ma^ 
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[ 12 ] 


) ^ M4 


(12)(13)(14)(23)(24)(34) 


(3 4) 

[ 12 ]- 


amplitude (96). 


With the identification (98), modulo the numerical coefficient, this matches the first term in the gravity 


For the second term in (96) that is obtained by a simple evaluation of a 4-vertex on shell, one 


should not expect the pattern (97) to continue. Indeed, the reason for si 2 in (97) is to cancel an 


extra internal line propagator when combining two YM amplitudes. In the case when an amplitude is 
obtained directly from a 4-vertex, with no internal propagator to be cancelled, one should expect the 
gravity amplitude to be simply the square of the YM one, as it happens for example for the 3-point 
amplitudes. And indeed, we see that the gravity all plus amplitude is a sum of three terms, each being 
a square of a colour-ordered gauge theory amplitude. The sum is necessary to produce a completely 
permutation-invariant expression from the ordered gauge theory amplitudes. 

So, overall we see that the gravity as a square of gauge theory pattern is at play for all the ampli¬ 
tudes. There are however subtleties about numerical coefficients still to be understood. Presumably, 
to reproduce the gravity amplitudes together with the correct numerical coefficient one should subtract 
the unwanted modes arising by tensoring the two spin one states; see m for a related discussion. It 
is also necessary to understand what replaces the usual double copy rules for the all plus amplitudes, 
with their two different types of terms. We leave all this to future work. 
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7 The first look at renormalisation 


As we have already mentioned in the Introduction, in many respects our theories behave as the 
renormalisable Yang-Mills theory from which they originated. This refers to the property of our 
perturbation theory that guarantees that there is at most a single derivative on every internal Feynman 
diagram line. In this aection we would like to analyse consequences of this for the power-counting, 
as well as study divergences of some simple diagrams. We will see that there are divergences of a 
new type, not present in YM, as can be expected from a power-counting non-renormalisable theory. 
However, at least at one loop these divergences can be absorbed by a local field redefinition. 


7.1 Choice of the renormalisation scheme 

As it was clear from our motivation section, and from the analysis of the tree-level scattering am¬ 
plitudes, the “collapsing” property (55) plays a very useful role in simplifying the structure of the 


perturbation theory. Thus, we definitely want to keep this property operating also at the loop level. 

However, here a complication arises. The most convenient renormalisation scheme to use in the 
context of such a theory as Yang-Mills is the dimensional regularisation. This is however not directly 
applicable to theories with spinors in them, because of the difficulties in continuing the notion of a 
spinor to arbitrary (non-integer) dimension. To put this problem differently, our theories use self-dual 
objects, and this is clearly a 4-dimensional notion, not directly extendible to arbitrary dimension. 

We will circumvent these problems in the following way. First, we note that we only need to 
dimensionally continue the loop momenta in our diagram, as it is the divergent integrals over these 
loop momenta that need to be regularised. So, the only object from the Feynman rules spelled out 
above that needs to be defined in general dimension is the momentum Uaa' on an internal line of a 
diagram. In four dimensions this object is defined as kAA' = where are the soldering 

forms, is the usual 4-momentum vector, and the minus sign is due to our choice of conventions, 
which are those of section 6 of [lOj . 

Now, in d dimensions one can define the notion of momentum kf^ in the usual way. To define 
the notion of kAA' we therefore need to dimensionally continue the soldering forms There is a 

known way to do this, familiar from the literature on the dimensional regularisation of theories with 
fermions. Indeed, a similar problem arises e.g. in QED, where one needs to dimensionally continue 
the 7 -matrices. One does this in a such a way as to preserve the Clifford algebra defining relation. 
With the soldering forms being just the 2x2 building blocks of the y-matrices, one can adopt a similar 
strategy to 0^^,. A practical way to do this is described in Appendix B.2 of [22], with some changes 
in conventions. In particular, it is explained in this work that some of the soldering form identities do 
continue to hold even in d 7 ^ 4 dimensions. The particular identity that is of interest to us is 


git^A e,q^)be' ^ 


(99) 


where now is the d-dimensional metric, while the spinor indices A, A' continue to be 2-dimensional. 
In particular = 2 as before. Contracting this identity with kuky we have 






( 100 ) 


which is the identity we have used in proving (55), but which now holds for a d-dimensional momentum. 


To summarise, our discussion has established that the identity (100) can still be used even in 
dimensional regularisation. One can thus continue to use the Feynman rules in their spinor form 
derived above. One should just be careful not to use in manipulations with the loop momenta some 
identity of the soldering forms that only holds in four dimensions. In calculations below we shall only 


use the identity (100), which is also what is required for the property (55) to hold. 
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It is important to stress that in our scheme the gauge field remains a 4-dimensional object, i.e. 
a spinor. In other words, we do not write and continue ^4“ to d dimensions. There is 

no need to do this in our setting as all Lagrangians and the resulting Feynman rules directly operate 
with objects kAA' and A\j^,. We need to continue kAA' to d dimensions for reasons of regularising the 
divergent integrals, but there is no reason to do the same with The propagator can continue to 

be given by (45), with the square of the d-dimensional k^ in the denominator, but the product of the 


usual 2-dimensional spinor metrics in the numerator. This is essentially the same scheme that is used 
in the 2-component version of QED in [22], with understanding that the field A°^j^i is now treated as 
a spinor, not tensor object. With these choices there is no difficulty in treating the considered here 
“self-dual” objects also in dimensional regularisation. 


7.2 There are never two dots on any internal line 


Having established that our identity (55) also holds in dimensionally regularised loop integrals, we 


would like to analyse consequences of this for the power counting. The identity (55) says that when 


two derivatives happen to act on the same internal line, this collapses the line. Importantly, this is 
as general as we stated it, because the derivatives of the gauge field only appear in the Lagrangian in 


the combination (10). In other words, they only appear as dots, and the property (55) applies. This 


means that there can never be more than one derivative on any internal line of any Feynman diagram 
that we need to analyse. We have already used this fact in our analysis of the large z behaviour 
of amplitudes under the BCFW shifts. Here we will analyse consequences of this fact for the power 
counting. 


7.3 Power-counting in effective field theory 

To understand the power counting, let us first recall the power-counting in a “typical” effective field 
theory. In this case we have Weinberg’s formula [23] 

( 101 ) 


D = 2L + 2 + — 2)mfc, 


where D is the superficial degree of divergence, L is the number of loops and is the number of 
vertices in the diagram with k derivatives in them. In case of many effective field theories, such as 
e.g. the chiral perturbation theory or gravity, there can only be an even number of derivatives at 


each vertex, and this number is greater or equal than two. Then the last term in (101) is a sum 


of non-negative terms. A particular case is that of perturbative gravity, where all vertices have two 


derivatives in them. In this case Hgr = 2L -|- 2. The relation (101) means that for each given D there 


is only a finite number of diagrams to consider. The relation (101) is easily obtained from the two 
obvious relations 


D = 4.L-2I + '^kmJ^, L = /-^m^-hl, 


( 102 ) 


where I is the number of internal lines. 


7.4 Power-counting in gauge theory 

In the case of gauge theories, such as theories of type f{F), with F being the (full or self-dual) 


curvature, the formula (101) still applies, but one also has vertices with one or no derivatives in them. 


Then the last term in (101) for /c = 0,1 contributes negatively. In such a situation this formula is not 


very useful. And indeed, in the case of usual YM theory, when we have just the cubic single derivative 
vertex and quartic no-derivative vertex we have an additional relation E + 21 = 3 V 3 -|- 414, which 
implies Dym = 4 — E typical of a renormalisable theory. 
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So, let us write down the power-counting relations for a gauge theory of the type f{F), first without 
taking into account the collapsing property (55). It is easy to see that in a theory of such type, at a 
vertex of valency n there can be as many as n derivatives, as well as n — 2, re — 4,..., i.e. the number 
of derivatives decreases by two. An example of this is the usual 1-derivative YM cubic vertex and the 
3-derivative new vertex. Thus, we label the number of re-valent vertices with k derivatives in them by 
We then have 


oo [^./S] 

D = AL-2I + J2J2(^- 

n=3 i=0 

oo [’^/2] 

i = '-EE + 1, (103) 

n=3 2=0 

oo [^/2] 

n=3 i=0 

where we indicated the summation limits explicitly. From the first and the last of these relations we 
get 


oo [>t./2] 

AL + E-D = Y,Y. (20C' > 0 , 

n=3 2=0 


(104) 


and thus the most divergent diagrams are those where all vertices have the maximal number of 
derivatives possible. In other words, we have 


D <4L + E, 


(105) 


with the equality being realised when all vertices have the maximal number of derivatives possible. 


It also should be noted that the maximal possible degree according to (105) is actually an overcount 
as far as the divergence is concerned. Indeed, to get the maximal possible degree of divergence we 
take all vertices to have the number of derivatives equal to their valency. It is then clear that precisely 
E of the derivatives are located on the external lines, and do not contribute to the divergence. Thus, 
the formula Dmax = ‘iL + E should actually be interpreted as a single quartic divergence per every 
loop integration. 

The degree of divergence that increases with the number of loops as 4L is clearly much worse than 
the divergence increasing as 2L, which is the case in e.g. GR. Below we will see that the collapsing 
property implies that the degree of divergence for our theories is in between these two cases. 

7.5 Power counting for deformations of YM 

We now take into account the property that when two derivatives act on a single internal line, the 


line collapses. This property implies that the divergence is never as bad as in (105). Indeed, to get 


an equality in (105) one needs to cancel all propagators by numerators, which precisely means two 


derivatives on each internal line. 

Let us write down relations for the power counting in the case when the maximal number of 
derivatives on each line is one. Let us denote by Ii the number of internal lines with a single derivative, 
and Jo those with no derivatives. We then have 


D — AL — Ii — 2 Iq + El , 

where Ei is the number of derivatives on the external lines. 


( 106 ) 
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Let us now compute the largest possible degree of divergence given the number of loops and the 
number of external lines. It is clear from (106) that the largest degree of divergence is achieved when 
all internal lines carry one derivative. So, let us assume that I = Ii. From the second relation in (103) 
we see that we minimise the number of internal lines at fixed L by having as small number of vertices 
as possible. Since loop diagrams with just one vertex vanish in dimensional regularisation, we need at 
least two vertices, which means that 


/min = L + l, (107) 

and 

L'max = 3L-l + F;i. (108) 

It remains to compute the maximal possible number of the derivatives on the external lines. It may 
seem at first that this number is equal to the number of external lines, but this is not true. Indeed, 
consider for example the diagram 



(109) 


It satisfies all the requirements assumed above: minimal number of vertices and thus maximal number 
of internal lines, each internal line carries single derivative. However, there is no quartic vertex with 
single derivative that could be used in this diagram, and so we learn that the maximal value of Ei 
cannot always be equal to E. The possible diagram in this case is 



( 110 ) 


and has the degree of divergence H = 2- 4 — 3 + 1 = 6. 

Simple analysis then shows that in the situation of minimal I the maximal number of dots on the 
external lines is E when L is odd, and E — 1 when L is even. Thus, we get 


f 3L — 1 + E, L odd 

\ 3L — 2 + E, L even 


( 111 ) 


In particular, this formula implies that Umax — E is always even. 

The formula (111) should be compared to (105) and (101). It is clear that the divergences in our 
case are better than (105), because the degree of divergence increases only as 3L and not 4L as would 
be the case in general theories of type f{F). On the other hand, we get a faster growth than e.g. in 
the case of perturbative gravity with D = 2L + 2. This is because there can be an arbitrary number of 
derivatives in a vertex in our case, as compared to just two in the case of the perturbative expansion 
of the Einstein-Hilbert action. 
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7.6 1-loop self-energy diagrams 


In order to understand if our class of theories has any hopes of being renormalisable, let us consider 
the simplest possible divergences arising at one loop. Since we already know how the divergences of 
YM theory manifest themselves, we only need to consider the modifications due to new interactions. 

Thus, let us consider the new self-energy diagrams that appear when the theory is deformed. We 
use dimensional regularisation, with the only subtlety being that the dimensionally continued objects 


kAA' satisfy only some of their properties in four dimensions. Above we have established that (100) 
can continue to be used. 

We do not consider tadpole-type diagrams that can be built from quartic vertices because these 
do not have any divergences in dimensional regularisation used here. This also applies to diagrams 
built from two cubic vertices where there is a pair of dots on one of the (or both) lines, because these 
diagrams collapse to tadpole-type ones. Thus, in particular there is no divergent diagram built from 


two vertices (50), because both internal lines of this diagram collapse. 


The only diagram whose divergence is non-zero in dimensional regularisation is one built from the 
new cubic and the YM cubic vertices, and where the dot of the YM vertex is on the external line 



( 112 ) 


It is clear that this diagram diverges quadratically. Since there are already two derivatives present on 
the external lines, it seems that the divergent part of this diagram requires a counterterm of type that 
is not present in the original Lagrangian. This is true, but just to an extent, because this counterterm 
corresponds to a simple field redefinition, as we shall now see. 

Using the Feynman rules, we see that to evaluate the above diagram we need to perform the Wick 
contraction in the following expression 

- k)A{i - k)rc^ (ns) 

where 1/2 is the symmetry factor. The Wick contraction gives 


We now replace I = q + xk, l — k = q + {x — l)k, and keep only the quadratic and zeroth order terms in 
the integration variable q. We also note that the primed indices of I and I — k are getting contracted. 
So, we will also need qsA'Qc^ = {^/‘^)^cbq‘^, which continues to hold even in d dimensions, as we have 


established in (100). With these replacements, the complicated structure appearing in (114) becomes 


^(g^ -k x{x - l)k‘^){4eM^€BN + es^eMN)- 


(115) 


The last term here is irrelevant as it contracts into the quantity that is AB and MN symmetric. The 
divergent part of the integral with (f' is two times the divergent part of the integral without q^ in the 
numerator, thus overall, after performing the x integration we get for this diagram 


1 \ag^C2 

(47r)2 


{kA{k)))\j,{kA{-k))) 


aAB 


(116) 
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This means that this divergence is removed by the foiiowing iocai heid redefinition 


A ^ ^ + 


1 ag^C2 k‘^ . 
e (47r)2 


appiied to the two point function —{i/2){kA{k)))‘\^{kA{—k)))°‘^^. 


( 117 ) 


7.7 General case 


Let us now try to generaiise the previous example. We know that at the one loop level, the maximal 
possible degree of divergence of a diagram is Umax = 2 + E. Since we already have in our Lagrangian 
Li-valent vertices with E derivatives in them, we see that there is a hope to renormalise the theory at 
one loop by removing the divergences of the maximal degree Umax = 2 + E using the field redefinition 
(117). Indeed, by applying this held redehnition to all the vertices of the original Lagrangian one may 
hope to remove all the maximal degree divergences. What remains are then the divergences of degree 
E, which are just the logarithmic divergences contributing to the multiplicative renormalisation of all 
the couplings. Below we will see how to realise this program at one loop for the renormalisation of 
the coupling constant a. 

We did not work out any details for the higher loop case, but if renormalisation is to work, it has 
to proceed by the local held redehnition of the type 


^ + A, (118) 

n=l ^ 2 

where Cn are divergent coefficients, followed by multiplicative renormalisation of all the couplings. 
One may also need non-linearities in the gauge held appearing in the held redehnition. Little more 
can be said without doing explicit calculations at higher loops, and we leave this to future work. 


8 /3-function calculation 

In dimensional regularisation (unlike in any Wilsonian scheme) only dimensionless couplings run. In 
our case all dimensionful couplings are expressed as multiples of powers of the same dimensionful 
parameter M. In dimensional regularisation this scale cannot run, with only the dimensionless cou¬ 
plings sitting in front of inverse powers of M receiving (inhnite) dimensionless corrections. This will 
be explicitly visible in the computations we perform, in the sense that the inverse powers of M present 
in the vertices always just go outside of the divergent integrals to be computed, and serve as the 
bookkeeping device for which term in the Lagrangian the divergence in question contributes to. To 
put it concisely, the scale M does not run in dim reg, and only the dimensionless couplings do. 

This feature of the scheme we use is quite important, for it implies that the presence of the new 
couplings does not change the running of the YM coupling g^. This decoupling of the flow of g^ from 
the flow for the other couplings does not occur in Wilsonian schemes. We find this decoupling nice 
and highly desirable, which serves as one more motivation to prefer the dimensional regnlarisation to 
any other possible scheme. 

Our strategy for computing the /3-function for the new coupling a is as follows. As we have seen 
previonsly, the presence of the new vertex does not change the multiplicative renormalisation of the 
gauge field from what it is in pure YM. Thus, we do not need to recompute the self-energy diagrams, 
as the result for the logarithmic divergence appearing in this case is known. We only need to consider 
the triangle diagrams. 

We then note that the presence of the new vertex cannot affect the logarithmic divergences con¬ 
tributing to the renormalisation of the usual YM vertex. This follows from the fact that the new vertex 
comes with a dimensionful prefactor of 1/M^. Thus, whatever the divergences that the diagrams with 
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the new vertex produce, they require counterterms of different mass dimension. Thus, there cannot 
be contributions to the running of the YM coupling g from the presence of the new interactions. 

The new interaction can thus only contribute to its own strength running. In order to find this in 
as simple way as possible, we evaluate the divergences on shell, for a configuration of three positive 
helicity gluons, in the same fashion as the new vertex was evaluated on-shell (70) in our discussion of 
the scattering amplitudes. We remind the reader that the three external momenta need to be continued 
to complex values to obtain a non-zero result. This on-shell evaluation simplifies the calculation and 
allows us to extract the divergent contributions to the new vertex rather easily. 


8.1 3-point triangle diagrams 


Let us consider possible triangle contributions to the vertex correction. We only need to analyse the 
effect of the new interactions, assuming the usual YM story as known. We note that diagrams with 
three new vertices completely collapse and are thus zero in dimensional regularisation. Diagrams with 
two new vertices become bubble diagrams proportional to the square of one of the external momenta, 
and hence are zero on-shell. Therefore we need only to consider the diagrams with one new cubic 
vertex. In this case there are three possible ways for each of the old vertex to place its dot, and thus 
in total nine diagrams to consider. In one of these diagrams the dots of both old vertices are placed 
on the internal lines already containing the dots of the new vertex. 


This diagram collapses and is 
those where no line collapses 



(119) 


regularisation. The remaining diagrams are 



( 120 ) 


as well as diagrams where a single line collapses 




( 121 ) 


plus an additional diagram 
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In this last diagram the bottom line collapses, and the result is a diagram proportional to k'^ and thus 
zero on shell. It remains to consider diagrams (120) and (121). 

Let us first consider (121). We need to compare the resulting collapsed diagrams to similar diagrams 
that appear if one pairs the old cubic and the new quartic vertex, see below. To evaluate the resulting 
collapsed diagrams we use the rule (54) which can be read as saying that whenever the two dots appear 
on the same internal line, the line collapses in the sense of being replaced by the identity matrix in the 
space of symmetric rank 2 spinors, times the numerical factor of 2i, which is what the numerical factor 
in (54) becomes if one takes into account an additional sign from the fact that one of the momenta is 
outgoing and the other incoming. Carefully comparing the corresponding Feynman rules we find 



and 





( 122 ) 


(123) 


as well as similar relations for the other collapsing line. 

We can now see that the diagrams in (121) precisely cancel similar diagrams from the old cubic 
new quartic vertex diagrams of the next section. Indeed, the symmetry factor in the triangle diagram 
case is 1/2. When the new vertex is on the line with the external momentum ki we get four collapsing 
diagrams. In two of them the momentum on the old cubic vertex of the effective bubble diagram is k^ 
and in the other two it is k 2 - When we take a sum over where the new cubic vertex can get attached we 
get overall four effective old cubic new quartic type diagrams with the momentum ki attached to the 
old cubic vertex. They all come with the symmetry factor of 1/2, and in two of them there is one dot 
inside the loop, and in the other two there are two dots. Overall we get precisely the same diagrams 
as appear in the last two terms in (130), with the same numerical factors, but with an opposite sign 
so these diagrams cancel with each other. 


The only triangles remaining are then the two types of diagrams in (120). The first one is 



{ag/2M^)f^^^(kiA^{ki)'^ (^{l - ki)A%l - h)) 

gfdef (^k2A\k2)) - k2)MM’Af{l)N^' (124) 

gfP<i(k^A^{k:i)Y'^AP{-l + k^)pp,A^{l + k2)Q^' 
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Wick contracting gives 


2M2 



r 


ade 


kiA^{ki))^^(k2A\k2)) 


MN 



PQ 


N (■ 


2{l — f"' P' ^MQ EM'' 


W' 


iZ2i(/ — A:i)2i(i; + k2)'^ 


Both the momenta in the numerator can be replaced by the loop integration momentum q = I — xki + 
yk 2 because this is the only divergent contribution. After this replacement is done the complicated 
product of e’s and the l,l — ki momenta reduces to 


2(?^eArBep"^eQM, 


plus a term proportional to eb^ which vanishes when contracted to the rest of the expression. The 
divergent part of this diagram is thus 


1 ag^C2 

2 (47r)2M2 



The other triangle diagram in (120) that we need to compute is 


A 


C 


(125) 



{ag/M‘^)f^^^(kiA^{kx)^ ^[-IA\-1)^ ^({I - ki)A%1 - ki)^ ^ 

gfdef _ k2)A<^{-l - k2))^^ A-{l)MN'Af{k2)N^' (126) 

gfP<l(^k3A^(k3)Y'^AP{-l + ki)pp:A1{l + k2)Q^' 


We do not include the symmetry factor of 1/2 because there are two diagrams with the dot being 
applied to either the /cs or k 2 legs. We shall see that on-shell the result is permutation symmetric, so 
adding the two diagrams cancels the symmetry factor. Performing the contractions gives 

- ^ rI’-[kiA'‘(ki)) {k:,A’-{k:^f'^ (127) 

- ki)fcc'P^p/2{l + k2)^^^ U2i(j — ki)^i{l + A:2)2 ’ 


where the minus in the brackets comes from the product of three structure constants and one extra 
overall minus comes from the fact that there are three propagators, each carrying a minus sign. We 
now need to open the symmetrisations, replacing the momenta with the loop integration momentum 
appropriately shifted, i.e. I = q + xki — yk 2 ,1 — ki = q + {x — l)ki — yk 2 ,l + k 2 = q + xki -|- (1 — y)k 2 - 
We then keep terms quadratic in q and replace q^qi, with {1/A)q^. Assuming the momenta ki,k 2 are 
on-shell, Iaa' = ^A^AA^AA' = ‘^^A^AA arid replacing the connections with the polarisation spinors 


using (69) gives the following result 


4(1 - 3y)(l 2) (2 3)(3 1) + 2(1 + 3x)(l 2)(3 


times a factor of (—i)^. Using the momentum conservation we can replace [pl]/[p2] = —(3 2)/(31). 
Thus, after the integration over the Feynman parameters we get 
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Thus, finally the divergent part of this diagram is 


1 ag^C2 2 
4 (47r)2M2 e 



( 128 ) 


8.2 3-point bubble diagrams: New cubic-new quartic/New cubic-old quartic 

As for bubble diagrams, hrst we note that any dots collapsing occurring in a bubble diagram can only 
lead to tadpole or simply to a point, both are known to be vanishing. In the presence of a new cubic 
vertex this leaves us only the following two diagrams. 




(129) 


However a straightforward calculation shows that both diagrams are proportional to momentum square 
and hence are zero on-shell. 

8.3 3-point bubble diagrams: Old cubic-new quartic 

There are three ways to attach the cubic vertex moving the dot and six ways to do the same with 
the quartic vertex. In total there are 18 diagrams. The symmetry factor for these diagrams is 1/4 
corresponding to flipping the two propagators, as well as flipping the external 23 legs. Among the 
resulting 18 diagrams 6 vanish because there is a double dot appearing on a propagator. Among the 
remaining ones there are just 5 different diagram topologies to consider 



For the last two diagrams it is understood that the four diagrams are actually 2 -|- 2 with two different 
arrangements for where the dot goes on the bottom legs. One more diagram is 



(131) 


However, this last diagram is proportional to kf and hence zero on shell. It will not be considered any 
further. As we have seen above, the last two diagrams in (130) are cancelled by what comes from the 
triangle diagrams. Hence, we only need to compute the first two diagrams. 
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The first one of them is very easy to compute 



{l/4)gr^'^(^kiA-{ki)yA’^{-l)AB'A‘^{l - ki)B 

o / \ / \ 


B' 


{-2i)^rP^r^(^k2AP{k2))^^{k3Ai{k5))^PA^{-l + ki)^pp,A\l)^^P' 


Here 1/4 is the symmetry factor. We now do the contractions of the placeholders for the internal 
gauge field lines. This results in the following integrand 


1 

4 ' 


{-2i)^{-C2rP<^)[kiA‘^{ki)y[k2AP{k2))^^{k3A<^{k3))/eBipe^'p'SA^^SB' 


1 


i/2i(/-A:i)2' 


The divergent part of the loop integral is then 
ag^C2 2 


-J^'^‘^{kiA^ 


B 




(47r)2M2 e" V ‘ ' 'J A 

The second diagram is a bit more difficult, but still rather straightforward 


(132) 



(1/4)5/“'’^((/ - ki)A%l - ki)^A\ki)AB'A^{-l)B^' 

{-2i)^rP<ir^(k2AP{k2))j^[kzA<i{k^))^PA^{-l + ki\pp,A\l)^)P' 


Doing the contractions gives 

.,3 




A' M)^B'P' 

epi ep 'e 

1 


h2i(^-A: i)2' 


To evaluate this we use the fact when I is replaced with I = q + xki only the ki term survives. Thus, 
after some algebra of the e’s, and using the on-shell property of Ai, we get 


(i)(x - 1)^ C 2 (^-l[k,A\k,))^^[k2AP{k2))^^{k,A<i{k,))^ 
with the divergent part evaluating to 


1 


N J iP[{l-ki)^'' 


1 ag^C2 Ij^abcl,^ Aa 


4 (47r)2M2 e 

8.4 /^-function 


B 




(133) 


Finally, now we have all of the pieces required for the /3-function calculation. Combining the results 


of (125), (128), (132) and (133) we get 


11 ^ 
2 “4“^^4/ ~~2 
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and thus the divergent part of the 3-point function, evaluated on-shell is equal to 


where we have multiplied the result by the 3! symmetry factor. In order to produce a finite result for 
the coupling, this divergence must be canceled by the following counterterm: 

Z„, = l+6|^t (135) 

Schematically, we can write the renormalised Lagrangian as 

— -|- Zgg^^^‘^A‘^{dA) + . (136) 

As remarked earlier, we absorb the self-energy bubble through a field redefinition, without writing the 
corresponding counterterm here. The mutiplicative renormalisation factor Za is therefore unchanged 
and given by its YM value 




10 ^ 1 

3 (47r)2 e 


(137) 


In addition we have the following relations between the renormalised and bare quantities 

ZaA? = Ag, ZagOtgg'^^'^ A? = {ag)oA^. 

From these two we get 

{ag)Q = agiAl’^ZagZ~^^‘^. 

Combining the results from equations ( |135 ) and (137) we have 

log(Z„,Z-='/ 2 ) = ^ 


and therefore 


(47r)2 e 


Q ^ ^(«ff)o ^ 1 dag e _ g^C 2 


91og(^) ag d\og{g) 2 (47r)2’ 

where we have used dg"^/d\og{gL) = —eg^ -|- .... We then arrive at a positive /3-function 

dag ag^C 2 

dlog{g) (47r)2 


(138) 

(139) 

(140) 

(141) 

(142) 


8.5 Interpretation 


The result (142) is our final result for the /3-function of the new coupling. It nicely does not contain 


any numerical factors at all. As could have been anticipated, this sign of the /3-function means that 
the new non-renormalisable interaction grows strong at high energies. 

We also note that the growth of ag is just logarithmic. Indeed, the coupling ag changes as 
( 5 r^)“^Ai. Qiven that the YM coupling changes with energy as 


g\A) = 


g^igo) 


1 -h {llC2/3{4:7ry)g‘^{fio)log{ijA/nl) 


(143) 
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we have 


{ag){ii) = {ag){go) (l + {IIC 2 m^'nf)g‘^{go) /^ll)) 


3/11 


(144) 


Thus, if we deform the YM theory slightly by adding a new coupling ag at some energy scale, the 
departure from YM grows stronger as the energy increases. The flow will no longer take one to the 
simple free UV hxed point of YM theory. This can be taken as a manifestation of the asymptotic 
safety scenario: If one wants to end up at a given fixed point (free one in this case), one should stay 
on the critical surface which in this case means ag equal to zero. As the new coupling becomes order 
unity, the formula (144) for its running can no longer be trusted. It is a very interesting question 


where the flow with non-zero {ag)o takes one in the UV. 


9 Discussion 


In this paper we have motivated and introduced a new family of power-counting non-renormalisable 
gauge theories in four space-time dimensions. The Lagrangian is constructed from powers of the 
self-dual part of the held strength. This means that in the Lorentzian signature the Lagrangian is 
complex, which makes physical interpretation difficult, but we nevertheless proceeded with analysis of 
the scattering amplitudes and the renormalisation. 

For the scattering amplitudes our analysis can be justihed by going to a (admittedly unphysical) 

setup of signature (-,++), where the Lagrangian becomes real. In this case one can still talk 

about null momenta, and thus set up the scattering theory with its on-shell scattering amplitudes. 
Our hrst main result is then that the scattering amplitudes for our family of theories continue to be 
determinable from lower order amplitudes by BCFW recursion relations. For the case of amplitudes 
with at least one positive and one negative helicity gluon one can use the standard BCFW recursion 
with only a pair of momenta getting shifted. For the case of the all plus amplitudes one needs to use 
the more involved all-momentum shift. Another complication in the all plus amplitude case is that 
only a part of the amplitude is determined by the recursion. There is also a part of the n-point all 
plus amplitude that is obtained simply by evaluating the n-valent vertex on-shell. This way a number 
of new couplings enter the recursion at every particle number n. Thus, as we have already said in 
the Introduction, our theories are constructible in the sense that one does not need a Lagrangian to 
compute the amplitudes. They can all be determined from simple building blocks using the rules of 
recursion. 

We note that the all momentum shift has already been used in the context of on-shell reconstruction 
of general held theory amplitudes in m- This work obtained a general condition for an applicability of 
the all-line shift to effective held theory amplitudes. Another related reference is |24] . which considered 
the non-linear sigma-model. In this work the all-line shift was used to set up a BCFW recursion for the 
Berends-Giele current, from which the amplitudes can then be obtained by setting the single off-shell 
leg on-shell. 

It is clear that the described particular version of constructibility of our theories (in the sense that 
the amplitudes are determined from lower amplitudes by a recursion) relies on the special property 


(55) of the associated perturbation theory. It is this property that guaranteed the 1/z behaviour of 
the amplitudes under the BCFW and Risager’s shifts. So, without this “collapsing” property many of 
the constructions of this paper are simply inapplicable. In particular, the usual 2-line BCFW shift is 


only available due to the property (55). However, the applicability of the Risager all-line shift is more 


general, as is explained in m- Some amplitudes in a general effective field theory for the gauge field 
are constructible by the all-line shift; see [20] . 

We have also seen that the double copy structure of gravity amplitudes continues to hold even 
for ’’deformed” theories. Thus, we have seen that the graviton scattering amplitudes for ’’deformed” 
General Relativity continue to be given by squares of ’’deformed” YM amplitudes. This observation 
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is not entirely new. In the context of parity-invariant gauge and gravity theories augmented by higher 
derivative operators it has been discussed in e.g. [20], and in much more details in m- 

The usual explanation of the double copy pattern for gravity is the open/closed string duality. 
For some non-renormalisable operators this explanation is still valid, see m- However, it appears 
unlikely that our chiral models with their infinite numbers of independent couplings can be realised 
in string theory. So, if the double copy pattern can be shown to hold in the context of theories with 
inhnite number of independent couplings, then the open/closed string duality cannot be what is at 
the root of this pattern. Then the only other possible justification of the gravity equals YM squared 
behaviour would be the fact that the amplitudes in both cases are constructed using the same recursion 
relation from the building blocks that themselves follow the double copy rule. This has already been 
identified as the explanation of the double copy pattern in the case of usual GR and YM in |25|; see 
also references therein. So, if the double copy pattern indeed extends to an infinite-parametric context 
as is suggested by our results, then one would be able to conclude that it is the constructible nature 
of the amplitudes that is at the heart of the corresponding double copy structure. 

As far as loop computations are concerned, we have seen that it is the property (55) that leads 
to an improved UV behaviour of our theories, in the sense that the degree of divergence grows slower 
than it would for a theory with as many as n possible derivatives in an n-valent vertex. Given that in 
our case there is at most a single derivative on each internal Feynman diagram line, the divergences 
of our theories are in a sense similar to those of the usual YM theory. This does not mean that there 
are no divergences of a new type. But we have seen that at least at the one loop level, the new 
divergences can be absorbed into a simple local field redefinition. After this is done, there remain only 
logarithmic divergences contributing to the multiplicative renormalisation of the coupling constants. 
We then explicitly determined the resulting renormalisation group running of the coupling ag, and 
found a positive /S-function (142). This is the second main result of this work. 

In the accompanying paper |3| we perform the complete one-loop analysis of our theories using the 
background field method. We will see that our expectation that at one loop the theory is renormalisable 
by a local held redehnition, as well as by multiplicative renormalisations, is realised. We will also 
explicitly compute the arising one loop beta-function. The result (142) is then conhrmed by a different 
calculation. 

Finally, let us discuss possible directions for future research. The most interesting question is 
whether the described here family of theories continues to be closed under renormalisation at higher 
loops. It is difficult to answer this without doing explicit calculations, and we leave this to the future. 

A related remark is that the open issues of reality of the Lagrangian have little to do with questions 
of renormalisability. Indeed, renormalisation is customarily studied in the Euclidean signature setup. 
In our case this would render the Lagrangian real, with the Euclidean path internal thus being of a non- 
oscillatory type. Thus, one can study divergences and renormalisation even without understanding 
issues related to the Lagrangian becoming complex for the Lorentzian signature. 

Another important question is whether one can make sense of the theory in the Lorentzian sig¬ 
nature, by finding a positive-dehnite inner product that would make the theory unitary. This is a 
difficult question, probably requiring a better yet understanding of the finite-dimensional examples of 
the PT-symmetric quantum mechanics. 


Note added: Supersymmetry 

As suggested by an anonymous referee of this paper, we add a note addressing the question of whether 
it is possible to supersymmetrise our class of theories (0. We will only consider the case of iV = 1 
supersymmetry. 

In four space-time dimensions, = 1 super-YM theory is a theory of a gauge-held as well as 
a Lie algebra valued Majorana spinor. In 2-component spinor notations, a Majorana spinor is the 
Dirac spinor composed of a 2-component unprimed spinor and its complex conjugate primed one. The 
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generator of supersymmetry is another Majorana spinor. 

Given that our approach is chiral, with one type of spinors playing a more important role than the 
other, and also that the requirement of Hermiticity of the Lagrangian is dropped, we can consider the 
following version of = 1 supersymmetric YM. The Lagrangian is given by 

c. ^+ Da.U% ( 145 ) 

The supersymmetry transformations are as follows 

^'^A — ^AB^^ ■ ( 146 ) 


Here is a (constant) SUSY generator. Note that in the above version of SUSY, only the unprimed 
Lie algebra valued spinor is transformed, with the primed one not changing under the supersymmetry 
transformation. Thus, here the unpr imed and primed '0^ spinors are not complex conjugates of 
each other, and the Lagrangian (145) is not Hermitian. Using = Dai^5A°'^'^ and integrating 


by parts in the first term, it is easy to see that the Lagrangian (145) is invariant under (146), modulo 


a surface term. The term cubic in the Lie algebra valued spinor that results from the variation of the 


covariant derivative in the second term in (145) vanishes due to the Grassmann nature of the spinors 

Let us now consider a more general Lagrangian of the type studied in this paper, still supplemented 
with a kinetic term for a Lie algebra valued spinor 


c = f{F%s)+r^'DA'^rA- 


(147) 


Let us consider the supersymmetry variation of the first term in this Lagrangian, with the transfor¬ 
mation rule for the gauge field unchanged from the one in (146). We get 


^fiFAB) = {f'TABDA'^i^' 


aA'B\ 


(148) 


where (/Oab is the first derivative of the function / with respect to its argument. It is then easy to 
see that to make ( 147| ) SUSY-invariant we just have to change the definition of the transformation 
rule for the spinor 


^i’A — if )ab^ 


B 


(149) 


Thus, the bottom line is that our theories are as supersymmetrisable as the usual Yang-Mills, with 
only the transformation rule for the spinor needing a modification. 
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